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ABSTRACT 

We have investigated the properties of the stellar mass function in 
the globular cluster NGC 6397 through the use of a large set of HST 
observations. The latter include existing WFPC 2 images in the V and 
I bands, obtained at ~ 4^5 and 10' radial distances, as well as a series of 
deep images in the J and H bands obtained with the NIC 2 and NIC 3 
cameras of the NICMOS instrument pointed, respectively, to regions lo- 
cated ~ 4^5 and ~ 3'2 from the center. These observations span the 
region from ~ 1 to ~ 3 times the cluster's half-light radius (r m — 3'), 
and have been subjected to the same, homogeneous data processing so 
as to guarantee that the ensuing results could be directly compared to 
one another. We have built color-magnitude diagrams that we use to 
measure the luminosity function of main sequence stars extending from 
just below the turn-off all the way down to the Hydrogen burning limit. 
All luminosity functions derived in this way show the same, consistent 
behavior in that they all increase with decreasing luminosity up to a peak 
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at Mj ~ 8.5 or M H ~ 7 and then drop precipitously well before pho- 
tometric incompleteness becomes significant. Within the observational 
uncertainties, at Mj ~ 12 or Mh — 10.5 (~ 0.09 M ) the luminosity 
functions are compatible with zero. The direct comparison of our NIC 2 
field with previous WFPC 2 observations of the same area shows that 
down to Mh — 11 there are no more faint, red stars than those al- 
ready detected by the WFPC 2, thus excluding a significant population 
of faint, low-mass stars at the bottom of the main sequence. By applying 
the best available mass-luminosity relation appropriate to the metallic- 
ity of NGC 6397 and consistent with our color-magnitude diagrams to 
both the optical and IR data, we obtain a mass function that shows a 
break in slope at ~ 0.3 M . No single exponent power-law distribution 
is compatible with these data, regardless of the value of the exponent. 
We find that a dynamical model of the cluster can simultaneously repro- 
duce the luminosity functions observed in the core, at ~ 3'2, 4'5, and 10' 
away from the center, as well as the surface brightness and velocity dis- 
persion profiles of red giant stars only if the model IMF rises as m - L6±0 - 2 
in the range 0.8 — 0.3 M and then drops as m - 2 ^- 1 below ~ 0.3 M . 
Adopting a more physical log-normal distribution for the IMF, all these 
data taken together imply a best fit distribution with m c ~ 0.3 and 
a ~ 1.8. 

Subject headings: stars: Hertzsprung-Russel (HR) and C-M diagrams - 
stars: luminosity function, mass function - Galaxy: globular clusters: 
individual: NGC 6397 
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1. Introduction 

Among the Galactic globular clusters (GC), 
NGC6397 has so far been the preferred tar- 
get in the quest for the lowest mass stars 
that could tell us much about the IMF of 
these systems and clarify whether or not they 
can hide a substantial fraction of their total 
mass in the form of very light objects. The 
close proximity of NGC 6397 (~ 2 kpc) cou- 
pled with its low metallicity ([Fe/H]~ —1.9; 
Djorgovski 1993) makes stars at the bottom 
of the main sequence (MS) relatively easy to 
observe with the powerful cameras on board 
the HST, and many have undertaken this 
challenge. Our original work on this subject 
(Paresce, De Marchi, & Romaniello 1995) re- 
vealed for the first time a marked deficiency 
of low-mass stars with respect to expectations 
based on a power-law mass function (MF) in- 
creasing all the way to the Hydrogen-burning 
limit (Fahlman et al. 1989). These findings 
were questioned by Cool, Piotto, & King 1996 
and Chabrier & Mera (1997), but have since 
been largely confirmed by independent mea- 
surements (Mould et al. 1996; King et al. 
1998) which all concur to suggest a MF that 
flattens out below ~ 0.3 M . We argue in this 
paper that below 0.3 M the MF drops all the 
way to the detection limit at the bottom of 
the stellar MS. 

Although the results mentioned above seem 
to be consistent with each other, and thus 
rather robust, there are at least two practi- 
cal issues that need to be carefully addressed. 
First of all, it is fair to wonder to which extent 
this deficiency of low-mass stars at the bot- 
tom of the MS is physical and genuine and 
how much of it should rather be attributed 
to some limitation inherent in the instrument 
or method of data analysis used. In fact, 
the observations carried out so far are con- 



fined to the optical domain (V and / bands), 
while low- and very-low mass stars such as 
those near the H-burning limit are expected 
to have a spectral energy distribution that 
peaks in the near-IR at ~ 1.2 /im (Allard & 
Hauschildt 1995). If a large population of 
very faint and redder stars were present in this 
cluster, they might escape detection with the 
WFPC2. Secondly, the accurate determina- 
tion of a luminosity function (LF) from these 
observations requires a reliable assessment of 
the contamination due to background unre- 
solved galaxies and to field stars which are 
rather numerous at the low galactic latitude 
of NGC 6397 (b u ~ -12). Attempts have 
been made to account and correct for this con- 
tamination using both color and proper mo- 
tion information, but both methods are sub- 
ject to some uncertainties. 

Finally, a third, more general question needs 
to be addressed concerning the real meaning 
of these localized MF, namely whether the 
observed flattening at low masses is a fea- 
ture inherent in the global MF (GMF) or, 
even, in the IMF of the cluster, or whether 
it simply results from the internal dynamical 
modification of an originally completely dif- 
ferent mass distribution. Although Richer et 
al. (1991) have proposed that, if measured 
near the half-mass radius, the present day MF 
should closely resemble the IMF, no direct 
confirmation of this hypothesis exists yet. 

In this paper, we address all three issues 
above. The first two are best attacked with 
the new NICMOS camera on board the HST, 
which we have used to investigate the stel- 
lar population of NGC 6397 precisely at the 
wavelengths where low mass stars are ex- 
pected to shine the most, and to sample a 
control field located well outside of the clus- 
ter tidal radius to provide an independent 
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estimate of field object contamination. We 
address the third issue by modelling the dy- 
namical state of NGC 6397 in Section 5, where 
we put together all the LF currently mea- 
sured for this cluster at various radial dis- 
tances and compare them with the predictions 
of a multi-mass model that takes internal dy- 
namical evolution (mass segregation) into ac- 
count. 

2. The Data 

The results of this paper are based on 
observations obtained by us with both the 
NIC 2 and NIC 3 cameras of the NICMOS in- 
strument on board the HST (MacKenty et 
al. 1997). NIC 2 was used on 1997, Sep 
30 to observe a region located ~ 4f6NW 
of the center of NGC 6397 (hereafter field 
F2) through the F110W and F160W broad 
band filters with a total exposure time of 
10,235 s and 4,027s, respectively. The field of 
view of the NIC 2 camera (~ 19" xl9") was 
placed over an area previously imaged with 
the WFPC2 (see Paresce et al. 1995). A field 
closer to the cluster center was observed on 
1998, June 17 with the NIC 3 camera (during 
the second "NIC 3 Campaign"), when the ~ 
50" x 50"FOV was pointed to RA=17:40:54, 
DEC=-53:42:57 (J2000) corresponding to a 
distance of 3!2 from the center (hereafter 
Field F 1) and images were taken for a total 
exposure duration of 2,303 s in F110W and of 
703 s in F160W. In units of the cluster's half- 
light radius (r hl = 175" Djorgovski 1993), our 
NIC 2 and NIC 3 fields are located, respec- 
tively, at ~ 1.6 ru and ~ 1.1 ru- We also 
used the NIC 3 camera to observe a compar- 
ison field (hereafter SKY), located 2°NE of 
the cluster at RA=17:47:52, DEC=-51:41:17 
(J2000), to estimate the contamination due 
to field stars and correct for it in a statisti- 



cal way. The filters and exposure times are 
identical to those used for the field F 1. 

The quality of all these data sets is excel- 
lent, with PSF full width at half maximum 
(FWHM) of ~ 0"13 for NIC 2 and of ~ 0"19 
for NIC 3 (achieved through the reposition- 
ing of HST's secondary mirror operated dur- 
ing the NIC 3 campaign). Images were sub- 
jected to the standard STSDAS calibration 
(CALNICA) to remove the instrumental sig- 
nature (bias and flat field correction) and, in 
the case of the NIC 2 images, multiple frames 
of the same area and taken through the same 
filter were combined to improve the statistics. 
The average per-pixel count-rate of the back- 
ground is of order 0.03 ± 0.015 count s _1 in 
F110W and 0.02 ± 0.008 count s" 1 in F160W 
for the NIC 2 data, and 0.11 ± 0.01 count s" 1 
in F110W and 0.12 ± 0.01 count s" 1 in F160W 
for NIC 3 images. These values, and par- 
ticularly the standard deviation of the back- 
ground, are larger than those predicted with 
the NICMOS simulator, but still within a fac- 
tor of two of them. 

The automated star detection routine daophot. daofind 
was applied to the data, with a detection 
threshold conservatively set at 5 a above the 
local average background level (the threshold 
was set to 10<r for the NIC 3 cluster field as 
it is more crowded than the others, and the 
severe PSF under-sampling of NIC 3 makes 
it very difficult to detect and measure faint 
stars near brighter objects). We carefully 
examined by eye each individual object de- 
tected by daofind and discarded saturated 
stars, a number of features (PSF tendrils, 
noise spikes, etc.) that daofind had inter- 
preted as stars, as well few extended 
objects (with FWHM larger than twice that 
typical of point sources). The number of well 
defined objects detected in this way amounts 
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to 517 in the F 1 field, 58 in the F 2 field, and 
119 in the SKY comparison field. We should 
notice here that we have used our knowledge 
of the quantum efficiency and filter transmis- 
sion of NICMOS in the F110W and F160W 
filters (see Pulone et al. 1998) to select expo- 
sure times in these bands such that a typical 
M-type dwarf in NGC 6397 would be detected 
with a similar SNR in both bands. And in- 
deed, with an exposure time about ~ 3 times 
longer in F110W than in F160W, all objects 
detected in one band are normally also visible 
in the other. 

Since crowding is not too severe in any of 
our images, stellar fluxes were measured using 
the standard digiphot.apphot IRAF aperture 
photometry routine, following the prescrip- 
tion of the "core aperture photometry" tech- 
nique described in De Marchi et al. (1993). 
Instrumental magnitudes were then calibrated 
and converted into the HST magnitude sys- 
tem (STMAG) following the relation: 

m ST = -2.5 log (^j (1) 

where c is the number of counts measured 
for each star, U the inverse sensitivity of the 
instrumental setup (camera + filters), e the 
encircled energy (i.e. the aperture correc- 
tion), and t the total exposure time. The 
internal accuracy of our photometry ranges 
from < 0.05 mag at m 160 ~ 18 to ~ 0.5 mag 
at m 160 — 26, although, particularly for the 
F110W observations taken with NIC 3, the 
highly variable pixel response function of the 
camera can introduce statistical uncertainty 
of up to ~ 10% (Storrs 1998). Due to the 
current uncertainty affecting the photomet- 
ric zero points of the NIC 3 camera (Nota et 
al. 1998), our absolute photometry is accu- 
rate to within the 5% — 10% level. In order 



to make it easier to directly compare our re- 
sults with those obtained from the ground, we 
translate our measurements from the STMAG 
system into the VEGAMAG photometric sys- 
tem of the HST, defined as one in which the 
magnitude of Vega would be in all bands 
and which is more similar to the classical 
Johnson-Cousin ground-based system. We 
adopt the VEGAMAG system in this paper 
and do so by subtracting the zero-point con- 
stants of 2.3 mag and 3.7 mag from the values 
of the STMAG magnitudes in the F110W and 
F160W bands, respectively. In the following, 
we refer to these magnitudes as to J and H, 
respectively. 

With the magnitudes measured in this way 
we have produced the color-magnitude dia- 
grams (CMD) shown in Figure 1 (panels a - 
c). The MS of cluster stars is well defined 
in the Fl field (NIC 3, Figure la), extending 
from H ~ 14.5 to H ~ 22 where it broad- 
ens due to increasing photometric errors. The 
MS of the stars in the F 2 field (Figure 1 b) 
is sparser and less well defined, but it occu- 
pies the same region in the CMD as that of 
the F 1 field. The latter, however, appears 
consistently broader than the former at any 
magnitude. This effect is likely to be the re- 
sult of the non uniformity of the sensitivity 
across the surface of the pixels of the NIC 3 
camera (Storrs 1998), coupled with the severe 
under-sampling of the PSF. The expected en- 
suing uncertainty in the photometry is larger 
for the F110W data (of order 11 %) than for 
the F160W filter (~ 6 %), and is in agreement 
with the observed MS broadening. 

The lines in Figure 1 a and 1 b represent 
the theoretical MS of a lOGyr old popula- 
tion with the metallicity and distance appro- 
priate to NGC 6397 as obtained using the 
models of Baraffe et al. (1997). We have 
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assumed the metal content [Fe/H] = —1.9 
and extinction coefficient E(B — V) — 0.18 
as given by Djorgovski (1993), but have let 
the distance take on two different values, i.e. 
(m - M) = 11.7 from Djorgovski (1993; 
solid line in Figure 1) and (m — M) Q = 12.2 
from Reid & Gizis (1998; dashed line). The 
theoretical sequences reproduce very well our 
observations, and, particularly, the change of 
slope occurring at H ~ 17.5 and due to the 
change of the main source of opacity in the 
atmosphere of the stars from H~ to H2. Be- 
cause of the uncertainty in the zero point of 
our photometry and of the broadening of the 
MS due to the intra-pixel sensitivity variation 
of the NIC 3 camera we cannot decide which 
of the two distance scales is more appropriate, 
in that they are both consistent with our ob- 
servations. The mass values printed along the 
right-hand axes of Figure 1 are derived from 
the M-L relation of Baraffe et al. (1997) for a 
distance modulus of (m — M)o = 12.2 (Reid 
& Gizis 1998). 

Comparing the CMD of the F 1 field with 
that of the SKY field (Figure 1 c) immedi- 
ately reveals that field stars (mostly disk ob- 
jects, owing to the low Galactic latitude of 
NGC6397, 6 n = — 12) occupy exactly the 
same location in the diagram as the cluster 
objects. It is thus not possible to assign clus- 
ter membership to a star on the basis of its 
color alone. This is precisely why we have se- 
cured observations of the SKY field so as to 
be able to account in a statistical way for the 
contamination due to field stars and affecting 
our LF (see Section 3). 

A robust determination of the LF also re- 
quires proper correction for the photometric 
incompleteness due to crowding. We have 
thus carried out several standard artificial 
star tests by adding artificial stars with mag- 



nitude values covering the whole range spanned 
by our observations. For fields F2 and SKY 
we have generated 10 artificial images per fil- 
ter adding each time respectively 8 and 20 
stars, and have found that more than 90 % of 
the artificially added objects are recovered at 
any magnitude level. The limited crowding 
and our conservative choice of a 5 a detection 
threshold combine in such a way that artificial 
stars are always retrieved unless they fall onto 
another object. In field Fl, however, crowd- 
ing is more severe. We ran 20 tests per filter 
adding 40 stars each time and found that 90 % 
of the artificially added objects are recovered 
at H ~ 19 but the completeness gradually 
drops to ~ 65 % at H ~ 22 as shown in Ta- 
ble 1. 

To derive the LF of MS stars in the clus- 
ter we have counted the number of objects as 
a function of the H magnitude in the CMD 
corresponding to fields F 1 and F 2. We have 
determined in this way the number of stars in 
each 0.5 mag bin (1-mag bins were used for 
F2), which is then multiplied by the incom- 
pleteness coefficient. This latter step needs 
only be performed for stars in field F 1, as 
the photometric incompleteness is negligible 
in fields F2 and SKY, where the Poisson 

Table 1: Table 1 - Photometric completeness 



in the F 1 field 


H magnitude 


completeness 


15 


100% 


16 


99% 


17 


97% 


18 


94% 


19 


90% 


20 


88% 


21 


75% 


22 


67% 
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statistics of the counting process dominates 
the final uncertainty on the LF. In order to 
account in a statistical way for field star con- 
tamination, we then subtracted from each bin 
the number of objects measured in the same 
magnitude bin of the CMD of field SKY (the 
area covered by SKY is the same as F 1 and 
seven times larger than F 2). The LF obtained 
in this way are shown in Figure 2 in a linear 
scale as a function of the observed H band 
magnitude. 

Although the LF in field F2 (bottom of 
Figure 2) spans a narrower magnitude range 
than that obtained with NIC 3 in F 1 (top 
histogram) because of saturation of the few 
bright stars in field F2, within the experi- 
mental errors the two LF agree very well with 
each other over the common magnitude inter- 
val, once properly scaled. Indeed, although 
the former LF is affected by large uncertainty, 
its shape is not different from that measured 
in F 1, with which it agrees well after hav- 
ing been scaled vertically by a factor of ~ 5 
to account for the the larger area covered by 
Fl. Both LF show a peak near H ~ 19 
followed by a clear drop all the way down 
to H ~ 22, where they become statistically 
consistent with zero. The error bars plotted 
in Figure 2 reflect the ±1 a uncertainty ob- 
tained by combining the Poisson statistics on 
the number counts with the statistical uncer- 
tainty associated with the background sub- 
traction. We note here, however, that we 
have associated a 2 a error bar to the bins 
in which the net number of stars is < 5, as in 
these cases the uncertainty can largely exceed 
the determined value and one might prefer to 
use a 2 a upper limit instead. Nevertheless, 
it is noteworthy that the number of stars be- 
comes compatible with being zero at magni- 
tude very close to that expected for the MS 



H-burning limit for stars with the metallicity 
of NGC6397 (Baraffe et al. 1997). 

3. Comparing the Optical to the Infra- 
Red Data 

NGC6397 is the first globular cluster for 
which a deep, reliable LF became available 
from WFPC 2 photometry of MS stars all the 
way to the bottom of the MS (Paresce et al. 
1995). So far, two other independent stud- 
ies exist with the WFPC 2 that confirm the 
deficiency of faint stars near the half-light ra- 
dius of this cluster discovered by Paresce et 
al. (Cool et al. 1996; King et al. 1998). A di- 
rect comparison of those observations, carried 
out at optical wavelengths, with our new IR 
data could clarify whether the WFPC 2 is effi- 
cient at detecting faint, red dwarfs or whether 
it has missed a large amount of the very low 
mass objects that could potentially account 
for a conspicuous fraction of the total mass of 
GC. 

The first step in this attempt is to com- 
pare deep optical (/ band) and IR (H band) 
images of the same area (field F2), as we do 
in Figure 3. It is immediately evident that all 
the objects detected with NICMOS in the H 
band with a ~ 170 min exposure (left panel) 
are also present in the ~ 90 min WFPC 2 im- 
age through the / band filter (F814W, right 
panel). The arrow indicates a star with My = 
13.8, M I = 12.4, Mj = 11.4, and M H = 10.7 
which, using the distance modulus of Reid & 
Gizis (1998) and the M-L relation of Baraffe 
et al. (1997) for [M/H]= -1.5, translates to a 
mass of <^ 0.09 M , namely that expected for 
an object close to the H-burning limit in this 
cluster (Baraffe et al. 1997). Since this ob- 
ject is clearly visible in all bands (V,I,J, and 
H), it is unlikely that the WFPC 2 has under- 
estimated the number of low-mass stars. In 
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principle, then, one should expect the two in- 
struments to give the same LF, but since the 
number of stars common to both NICMOS 
and WFPC 2 is small, statistical fluctuations 
might be important. 

A second, perhaps better way to reveal any 
significant difference between the results of 
NICMOS and the WFPC 2 data is to com- 
pare the IR and optical LF to one another. 
The wavelength difference, however, would re- 
quire a transformation from one photometric 
system to the other, which would inevitably 
introduce uncertainties that we want to avoid. 
To minimize possible systematic errors, we 
prefer to compare to one another the theo- 
retical MF that best fit our visible and IR 
data. 

In Figure 4 a we show the LF of field F 1 in 
units of absolute if -band magnitude Mr, as- 
suming a distance modulus of m — M = 12.2 
(Reid & Gizis 1998) and E(B - V) = 0.18 
(Djorgovski 1993). Figure 4b shows the LF in 
the I band (F814W) as measured by Paresce 
et al. (1995) with the WFPC 2, also converted 
to absolute magnitude units Mj. The LF in 
the I band does not cover the whole MS, in 
that it stops at the bright end ~ 2 mag be- 
fore reaching the turn-off (Mi ~ 4) due to 
saturation. 

Since it is often assumed for simplicity 
that, even over a wide mass range, the stellar 
MF is well represented by a power-law dis- 
tribution of the form dN/dm = m~ a (see, 
however, Scalo 1998 about the large range 
spanned by a), we initially adopt such a 
functional form for the theoretical MF.f] The 

3 Some authors prefer to express the power-law MF in 
logarithmic units, i.e. dN/d\ogm — mT x . The two 
exponents a and x are such that a = x + 1, with 
Salpeter's IMF taking on the value of a = 2.35 and 
x = 1.35, respectively. Throughout the rest of the 



dashed lines in Figure 4 show the theoretical 
LF obtained by multiplying a power-law dis- 
tribution by the derivative of the M-L rela- 
tion appropriate to NGC6397 (Baraffe et al. 
1997). This M-L relation is adopted here be- 
cause it is the only one presently available 
that takes into account all the stellar physics 
on the problem and that fits very well all 
available data on the CMD of globular clus- 
ters. The inspection of Figure 4 a and 4 b 
immediately reveals that a single-exponent 
power-law distribution is not a viable choice 
for the underlying MF at this position, re- 
gardless of the value of a. And indeed, while 
the bright portion of the LF up to the peak 
at M H — 6, Mi ~ 8 requires a ~ 1.4 (ac- 
ceptable fits are found for 0.8 < a < 1.6 
in the H band, and 1.2 < a < 1.6 in the I 
band), a much shallower exponent a ~ —0.4 
(-0.8 < a < -0.4 in H, -0.5 < a < -0.1 in 
I) is needed to fit the dropping portion of the 
LF beyond the peak. At the same time, Fig- 
ure 4 shows that the same type of MF applies 
to both the optical and IR data, i.e. one that 
flattens out and drops at ~ 0.3 M , close to 
the peak of the LF when adopting the M-L re- 
lation of Baraffe et al. with the distance mod- 
ulus and reddening values given above. In 
other words, the peak and turn over observed 
in all the LF measured so far for NGC6397 
near the half- mass radius (Paresce et al. 1995, 
De Marchi & Paresce 1997, King et al. 1998) 
cannot be attributed exclusively to the change 
of slope of the M-L relation at ~ 0.3 M Q . 

This result is fully consistent with the con- 
clusions of Paresce et al. (1995) and De 
Marchi & Paresce (1997), showing that they 
apply to the IR as well. We can, there- 
paper we will use the linear notation (a), except for 
Figure 6 which, as we explain, is best displayed in 
logarithmic units. 
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fore, exclude that near the half-light radius 
of NGC6397 the MF of the stellar population 
is represented by a single- exponent power-law 
distribution as advocated by Chabrier & Mera 
(1997) and Silvestri et al. (1998). Chabrier 
& Mera (1997) reach this conclusion by as- 
suming that the last four bins of Paresce et 
al. (1995) can be safely ignored as being 
seriously affected by incompleteness. This 
is now known not to be the case. Silvestri 
et al. (1998) reach this conclusion by us- 
ing a model that is inconsistent with the ob- 
served CMD of NGC6397 and other clus- 
ters. This is most likely due to their use of 
a grey-like atmosphere approximation which 
Chabrier & Baraffe (1997) have convincingly 
demonstrated to be wrong below ~ 0.6 M . 

4. Local Mass Function or Global Mass 
Function ? 

It is well known that, through the relax- 
ation process, stars in GC tend to reach en- 
ergy equipartition and, as a result, their mass 
spectrum can vary in time and space (Spitzer 
1987) so that, for instance, the MF measured 
in the core of a cluster in thermal equilib- 
rium should decrease with decreasing lumi- 
nosity. Paresce, De Marchi, & Jedrzejewski 
(1995), King, Sosin, & Cool (1995), and De 
Marchi & Paresce (1996) have indeed shown 
this to be precisely the case, in the mass 
range 0.8 — 0.5 M , respectively for 47Tuc, 
NGC6397, and M 15. At this juncture, it is 
fair to wonder whether the flattening below 
~ 0.3 M that we now observe in the MF of 
field F 1 and F 2 is an intrinsic feature of the 
IMF, or perhaps a characteristic signature of 
the global MF, or whether it simply is the re- 
sult of some local dynamical effects bearing 
no relation whatsoever with either the IMF 
or the global MF. Although conventional wis- 



dom and theoretical models suggest that away 
from both the cluster center and from its pe- 
riphery these effects should be rather modest 
(Richer et al. 1991; Vesperini & Heggie 1997), 
no unambiguous proof exists that this hy- 
pothesis is correct, nor does one exactly know 
where (if anywhere) the local MF should be 
most insensitive to dynamical modifications. 
Thus, since NGC6397 is probably in an ad- 
vanced evolutionary state and has probably 
already gone through the collapse phase, the 
two-body relaxation process that has strongly 
changed the properties of the stellar popula- 
tion in the cluster's core (King et al. 1995) 
might have worked its effects farther out in 
the periphery, possibly imparting non negli- 
gible modifications onto the MF in the re- 
gions that we observed. And indeed, a care- 
ful inspection of Figure 4 already reveals that, 
although similar, the MF that best fits the 
LF in field F 1 is not precisely the same that 
one requires to properly reproduce the LF ob- 
served in field F2 (see discussion above). 

In an attempt to clarify this issue and to 
put more stringent constraints on the shape 
of the cluster's global MF (and, possibly, on 
its IMF), we have computed the effect of the 
segregation mechanism in NGC6397 by us- 
ing a dynamical model of the cluster so that 
we could compare our observations with the 
predicted shapes of the LF at several radial 
positions. Besides the data obtained in field 
Fl (r ~ 3f2) with NIC 3 and in field F2 
(r ~ 4?5) with WFPC2 and NIC 2, we have 
included in this exercise the data taken with 
the WFPC2 by Mould et al. (1996) in a field 
located ~ 10' N of the cluster center. In or- 
der to ensure that the comparison of these 
latter data to our own observations is not 
hampered by possible differences in the way 
they were reduced, we have extracted from 
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the HST archive the images taken by Mould 
et al. (1996) and have subjected them to the 
same, homogeneous data processing that we 
have applied to our WFPC2 observations of 
NGC6397 (Paresce et al. 1995). 

Observations of Mould et al.'s field (here- 
after called F3) are available through both 
the F555W and F814W bands, with a total 
exposure time of 3, 600 s in each filter. This 
has allowed us to build a CMD in all ways 
similar to that of Mould et al. (see their Fig- 
ure 1), to which we have applied our "2.5 a- 
clipping" criterion to discriminate, in a sta- 
tistical way, the stars belonging to the cluster 
from field objects (see De Marchi & Paresce 
1995 for a detailed description of the method). 
The LF of MS stars determined in this way 
agrees reasonably well with that measured 
by Mould et al., in that both LF increase 
with decreasing luminosity up to a peak at 
m 8 i4 ~ 20.5 and then fall off. Theirs, how- 
ever, reveals a secondary peak at m 8 i4 — 22 
(Mi ~ 9.5) that we do not see and which is 
most likely the result of a different amount 
of correction for field star contamination. We 
would like to point out that the latter repre- 
sents the largest source of uncertainty in the 
determination of the LF in field F3, and that 
this effect is much more severe than that made 
by ignoring the photometric incompleteness 
as both Mould et al. and we have done. The 
resulting LF is shown in Figure 5 after having 
been converted into units of absolute magni- 
tude by adopting the distance modulus given 
above. 

5. Modelling Mass Segregation 

We have employed the multi-mass Michie- 
King models developed by Meylan (1987, 1988), 
and which were later adopted by Meylan & 
Mayor (1991) to study the dynamical proper- 



ties of NGC6397. Each model is character- 
ized by an IMF in the form of an exponen- 
tial dN/dm = m~ a , with a variable exponent 
a, and by four structural parameters describ- 
ing respectively the scale radius (r c ), the scale 
velocity (v s ), the central value of the dimen- 
sionless gravitational potential W Q , and the 
anisotropy radius (r a ). From the parameter 
space defined in this way, Meylan & Mayor 
(1991) have selected those models that simul- 
taneously fit both the observed surface bright- 
ness (SBP) and velocity dispersion (VDP) 
profiles of the cluster. The fit to the SBP and 
VDP, however, can only constrain r c , v s , W , 
and r a while still allowing the IMF to take on 
a variety of shapes. To break this degeneracy, 
we have complemented Meylan's code by im- 
posing the condition that the model MF agree 
with the observed LF, following precisely the 
approach that we used in our investigation of 
the globular cluster M4 (Pulone, De Marchi, 
& Paresce 1999). 

As we explain in that paper, Michie-King 
modeling only provides a "snapshot" of the 
current dynamical state of the cluster. It is, 
then, useful to define the global mass func- 
tion (GMF), i.e. the mass distribution of all 
cluster stars at present, as the MF that the 
cluster would have simply result of stel- 
lar evolution (i.e. ignoring any local modifi- 
cations induced by internal dynamics and/or 
the interaction with the Galactic tidal field). 
Clearly, in this case the IMF and GMF of MS 
(un-evolved) stars is the same. For practical 
purposes, the GMF has been divided into six- 
teen different mass classes, covering MS stars, 
white dwarfs, and heavy remnants, precisely 
as described in Pulone et al. (1999). 

The SBP and VDP used in our simulations 
are taken from the paper of Meylan & Mayor 
(1991), as G. Meylan kindly gave us the data 
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in electronic form. Since these authors pro- 
vide the mass, density, and structural param- 
eters of the 8 models that best reproduce the 
observed radial profiles (SBP and VDP), we 
have tested whether the predicted LF would 
also agree with our observations. Our exer- 
cise confirms what we have already shown in 
Figure 4: as long as a single value of the ex- 
ponent a is used for the IMF over the mass 
range 0.1 — 0.8 M , none of these 8 models 
gives local MF that, converted into LF using 
the M-L relation described above, agree with 
our data. In fact, if the exponent a is se- 
lected in such a way that the rising portion of 
the LF is well fitted (this requires a ~ 1.6), 
the model LF consistently deviate from the 
observations beyond the peak at Mi ~ 8.5, 
Mh — 6.5 because they fall off too slowly. 
Similarly, to fit the falling portion of the LF 
one would require a ~ —0.2, which would in 
turn overestimate the number of bright stars. 

On the other hand, if the model IMF is 
allowed to flatten out below ~ 0.3 M Q , both 
portions of the three observed LF can simul- 
taneously be fitted while, at the same time, 
the SBP and VDP are adequately reproduced. 
In particular, we find that an exponential 
IMF with index a = 2 ± 0.2 for evolved 
stars (m > 0.8 M ), a = 1.6 ± 0.2 in the 
range 0.8-0.3 M and which drops with a = 
—0.2 ±0.1 at lower masses would not only pro- 
duce a good fit to the LF observed in field F 1, 
F2, and F3 (Figure 5), but would also give 
reasonable values for the cluster's structural 
parameters: the concentration ratio (c ~ 2.6), 
the total cluster mass (~ 8.9 x 1O 4 M ), 
the fraction of heavy remnants (~ 33 % in 
mass, about half of which in the form of neu- 
tron stars of 1.2 M and the rest in the form 
of heavy white dwarfs of O.6M ), and the 
anisotropy radius (r a /r c ~ 20) are all in excel- 



lent agreement with those of the 8 best models 
of Meylan & Mayor (1991). 

King et al. (1998) notice that the some- 
what large value of the core radius that Mey- 
lan & Mayor (1991) obtain (and that we only 
marginally decrease here) does not allow them 
to accurately reproduce the SBP in the in- 
nermost regions of the cluster. On the other 
hand, Meylan & Mayor point out that the 
deviation of the observed profile from that 
of a typical Michie-King model in the cen- 
tral regions is most likely the result of ran- 
dom fluctuations due to a few bright stars, 
whose distribution is not uniform (Auriere, 
Lauzeral, & Ortolani 1990). Moreover, if the 
core of NGC 6397 has already been through 
its collapse phase, a classical King-type pro- 
file (King 1962) of a cluster in equilibrium cer- 
tainly does not provide the most appropriate 
representation of the SBP and the value of 
r c that one would deduce from it would be 
rather uncertain. 

To test whether our model also gives a 
proper account of the properties of the stel- 
lar population in the core, we have compared 
the MF that it predicts at ~ 7" from the cen- 
ter with that measured by King et al. (1995) 
with the FOC on board HST. Although our 
model seems to predict a somewhat slower de- 
cline than that reported by King et al. in the 
range 0.7 — 0.4 M (namely a ~ —2 vs. their 
a ~ —2.7), this discrepancy should not be 
taken too seriously in view of the limitations 
of our model and of the uncertainties that 
could affect the conversion of the observed LF 
from the F480LP band of the FOC into the 
Johnson's V band, and then from that into 
the F814W band of the WFPC 2 and, finally, 
into a MF by the use of a rather crude M- 
L relation. We can, therefore, conclude that 
the two MF are in fairly reasonable agreement 
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with each other, in that they both indicate 
an inverted function (i.e. stars are decreasing 
with decreasing mass) starting already at the 
turn-off. 

6. Discussion 

The result stemming from Figure 5 and 
from the discussion above is fully consistent 
with our findings in Section 3. Not only can 
we exclude that the local MF measured near 
the half-light radius come from a single expo- 
nent power-law, but we also prove here that 
not even the global MF from which they orig- 
inated through energy equipartition can be a 
function of that type. These findings carry 
the important implication that, if NGC6397 
has not experienced too strong an interaction 
with the Galactic tidal field or bulge (but 
see Gnedin & Ostriker 1997), also the IMF 
has to flatten out and, possibly, drop below 
~ 0.3 M . 

Since, however, Nature has always shown 
a predilection for smooth variations, it might 
seem rather unphysical that the IMF or GMF 
of MS stars can have such an abrupt change 
of slope at ~ 0.3 M as our work seems to 
imply. Clearly, our fit to the GMF is noth- 
ing but a mathematical over-simplification, 
dictated by the adoption of a power-law as 
representative of the distribution of stellar 
masses. As Scalo (1998) points out, this as- 
sumption might be correct for intermediate 
mass stars, but fails at the edges of the dis- 
tribution (very low and very high masses). 
Although we could have used more than two 
values of the exponent a to fit the change of 
slope, that would have appeared as an unwar- 
ranted sophistication whose physical meaning 
remains uncertain. Thus, since the most im- 
portant implication of Figures 4 and 5 is that 
the MF of NGC 6397 cannot be reproduced 



using a single exponent power-law distribu- 
tion, one wonders why choosing two or more 
power-laws would make any sense at all. 

There is, however, a large amount of the- 
oretical work suggesting that the star forma- 
tion process would preferentially give rise to 
a mass distribution that deviates from a pure 
power-law. Zinnecker (1984) and, more re- 
cently, Adams & Fatuzzo (1996) have pointed 
out that, if the number of parameters gov- 
erning the fragmentation process of proto- 
stellar clouds is sufficiently large (of order 5 or 
more), the central-limit theorem immediately 
implies that the IMF should take on a log- 
normal form. A similar IMF is obtained in- 
dependently by Elmegreen (1997, 1999) with 
a different theoretical approach. If we adopt 
the formalism of Adams & Fatuzzo (1996), 
the IMF is characterized by two parameters, 
namely the characteristic mass m c and the 
standard deviation a and takes on the form: 

h/OogmHA- ^y (2) 

where A is a normalization constant. The 
Miller & Scalo (1979) IMF (now superseded; 
Scalo 1998) would require m c ~ 0.095 M and 
a = 1.57, while the Scalo (1986) IMF (now 
also extinct; Scalo 1998) would imply both 
a larger value of m c and a smaller width a. 
In Figure 6 we show that, within the observa- 
tional uncertainties, the GMF of NGC 6397 is 
consistent with a log-normal distribution hav- 
ing m c ~ 0.3 M and a ~ 1.8. Interestingly, 
these values are very similar to those found 
for MS stars in the Pleiades (Bouvier et al. 
1998). 
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7. Summary and Conclusion 

The main results of this paper can be sum- 
marized as follows: 

• The properties of the stellar population 
of NGC6397 as determined with deep 
IR observations obtained with NICMOS 
on board the HST are fully consistent 
with those stemming from similarly deep 
WFPC2 optical data. In particular, 
down to the bottom of the cluster's MS, 
NICMOS does not show any more faint, 
red objects than those already detected 
by the WFPC2, thus confirming that 
there are no stars contributing signif- 
icantly to the low-mass budget of the 
cluster which can be detected in the IR 
and not in the / band. Thus, the un- 
availability of NICMOS will not hinder 
this investigation in the future. 

• A comparison of our deep IR CMD with 
the theoretical models of Baraffe et al. 
(1997) shows the agreement to be ex- 
cellent, thus proving that the latter are 
reliable also in the IR domain. 

• The IR LF of MS stars near the cluster's 
half-light radius, measured using a com- 
parison field located a few degrees away 
from the cluster center to correct for the 
contamination due to foreground and 
background objects, is fully consistent 
with those obtained with the WFPC2 
in the / band using different and inde- 
pendent techniques to separate cluster 
members from field stars, namely the 
color selection of Paresce et al. (1995) 
and the proper motion study of King et 
al. (1998). 

• Using the M-L relations of Baraffe et 
al. (1997), all three LF translate into 



MF in which the number of objects per 
unit mass increases exponentially as ~ 
m~ L5 in the range 0.7 — 0.3 M and then 
flattens out and drops. We show that 
a simple power-law distribution of the 
form m~ a cannot reproduce any of these 
LF if the exponent a is kept constant, 
regardless of its value. 

• The relaxation mechanism and the en- 
suing mass segregation process inside 
the cluster is studied under the assump- 
tion that a simple isothermal model of a 
cluster in equilibrium is a viable choice 
for NGC6397 and ignoring the effects 
of tidal stripping due to the Galaxy, as 
they do not seem to play any signifi- 
cant role in reshaping the MF of this 
cluster (Paresce & De Marchi 1999). 
The model can accurately and simulta- 
neously reproduce the LF observed in 
regions located ~ 3f2, 4f5, and 10' away 
from the center, the inverted LF in the 
core, as well as the SBP and VDP of 
red giant stars only if the IMF rises as 
m -i.6±o.2 j n ran g e o.8 — 0.3 M and 

then drops as rn°- 2±0A below ~ 0.3 M . 
The error bands associated with the 
quoted exponents assume that the M- 
L relations of Baraffe et al. (1997) are 
not affected by statistical uncertainties. 

• The IMF that we obtain is fully con- 
sistent with a log-normal distribution 
peaked at m c ~ 0.3 M and with stan- 
dard deviation a ~ 1.8. 

In this paper we show that a robust and 
reliable determination of the IMF of GC re- 
quires the effects of dynamical evolution to 
be taken into proper account. Although we 
have used the best LF currently available for 
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NGC6397, we are still forced to rely on un- 
certain brightness and velocity dispersion ra- 
dial profiles and on very simple Michie-King 
multi-mass quasi-equilibrium models to con- 
vert a locally determined present day mass 
function into a global MF. In fact, while the 
models that we have used are currently the 
state of the art, they work only under the 
strict condition that the cluster be fully re- 
laxed and in dynamical equilibrium. Episodes 
of core collapse and, especially, gravitational 
shocking due to repeated interactions with 
the bulge and disk of the Galaxy, however, 
may profoundly disrupt this distribution by 
ejecting low mass stars from the core and 
by compressing the tidal boundary in phase 
space at each encounter (see e.g. De Marchi 
et al. 1999). Until these effects can prop- 
erly be accounted for in future models of this 
sort, we should obtain the global cluster MF 
directly by observation through a systematic, 
precise, and complete census of the entire stel- 
lar population of GC spanning a wide range of 
possible evolutionary scenarios. Specifically, 
one needs to determine the exact shape of the 
cluster LF at many locations in the cluster 
from the main sequence turn-off down to at 
least the peak of the LF at ~ 0.3 M Q . These 
local LF and the resulting global MF can then 
be compared to the available cluster models 
to assess their validity under differing circum- 
stances. This project is within easy reach of 
a large telescope such as the VLT. 

We would like to thank Georges Meylan for 
providing us with the original code of his dy- 
namical model and Isabelle Baraffe and Gilles 
Chabrier for many useful discussions. 
We are indebted to Haldan Cohn, the referee 
of this work, for his invaluable comments that 
have considerably strengthened our paper. 
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Fig. 2. — IR luminosity functions as measured 
in field F 1 (top histogram) and F 2 (bottom 
histogram). Error bars reflect the ±1 a uncer- 
tainty of the Poisson statistics on the number 
counts and the statistical uncertainty associ- 
ated with the background subtraction. Note, 
however, that we have associated a 2 a error 
bar to the bins in which the net number of 
stars is < 5, as in these cases the uncertainty 
can largely exceed the determined value and 
one might prefer to use a 2 a upper limit in- 
stead. 



Fig. 1. — IR Color-magnitude diagrams as 
measured in field Fl (a), F2 (b), and SKY (c). 
The magnitudes are given in the VEGAMAG 16 
system. See the text for an explanation of the 
lines in panel (a) and (b). 



Fig. 3. — Comparison of the same area (field 
F2) as seen in a with NICMOS NIC 2 in 
the H band with a ~ 170 min exposure (left 
panel) and in a ~ 90 min duration WFPC2 
image through the / band filter (F814W, right 
panel). The arrow indicates a star close to the 
H-burning limit in this cluster (see text). 
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Fig. 4. — Comparison between the IR and op- 
tical LF: (a) LF of stars in field F 1 in units 
of absolute iJ-band magnitude Mr, and (b) 
the LF in the / band (F814W) as measured 
by Paresce et al. (1995) with the WFPC2, 
also converted to absolute magnitude units 
Mj. The dashed lines show the best fitting 
underlying MF, characterized by a = 1.2 and 
a = —0.6 respectively for the rising and drop- 
ping portion of the H band LF, and a = 1.4 
and a = —0.3 in the I band. Acceptable fits 
are obtained for a within about ±0.2 of the 
best fitting values, as explained in the text. 
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Fig. 5. — The solid lines show the LF pre- 
dicted by the dynamical model that best fits 
simultaneously the data at 3'.2 (triangles), 4f5 
(boxes), and 10' (diamonds) radial distance 
as well as the SBP and VDP (see text). The 
scale on the top axis refers only to the IR LF 
at r = 3'.2. The dashed line shows the LF as 
measured by King et al. (1998) in a different 
region of the cluster, yet still ~ 4'.5 away from 
the center. 
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Fig. 6.— Top panel: The GMF of NGC6397 
as predicted by our dynamical model (thick 
solid line) is consistent, within the error band 
marked by the thin lines, with a log-normal 
distribution having m c ~ 0.3 M Q and a ~ 1.8 
(dashed line). The error band corresponds to 
the uncertainty in the slopes of the two power- 
law segments used to represent the IMF in our 
dynamical model (see Sections 3 and 5) and, 
as such, it does not directly reflect the statis- 
tical fluctuations of each individual bin in our 
LF. Please note that, because the log-normal 
distribution is best represented in units of 
logarithmic mass, also the power-law curves 
in the figure are plotted in the same units, 

1. e. dN/dlogm = m~ x . The exponent a 
used throughout the text obeys the relation 
a = x + 1 as explained in Section 3. Bot- 
tom panel: Same GMF as above, but in lin- 
ear mass units, i.e. dN/dm. The dotted lines 
represent the MF of NGC 6397 as obtained by 
King et al. (1998) using the mass-luminosity 
relations of Baraffe et al. (1997), respectively 
for a distance modulus of 12.05 (boxes) and 
12.05 (triangles). See King et al. (1998) for 
the details. 




